Collinear ordering of easy-axis triangular lattice antiferromagnets 
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Antiferromagnetically coupled moments on the frustrated triangular lattice typically order in a 
coplanar state at low temperature. Here, we demonstrate that the presence of not-very-large easy 
axis single ion anisotropy leads to an interesting orientationally ordered collinear state in triangular 
lattice antiferromagnets with moments S > 3/2. This ordered state breaks the symmetry of n/3 
rotations about a lattice site, while leaving intact the translational symmetry of the lattice. 
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Introduction: Magnetic ions in insulating solids of- 
ten interact through 'exchange' couplings between near- 
est neighbours. In some cases, these exchange couplings 
compete with each other due to the geometry of the lat- 
tice. Such 'geometrically frustrated' magnets typically 
have a broad cooperative paramagnetic regime [l| for a 
range of temperatures below the exchange energy scale. 
The behaviour of the system in this regime is controlled 
by the manner in which the spins explore the multitude 
of classical ground states that arise from the geometric 
frustration. This gives way at very low temperature to a 
variety of novel phases [2| arising from quantum effects 
and sub-leading interactions. 

Antiferromagnets on the triangular lattice provide par- 
ticularly well-studied examples of magnetic frustration. 
When spins on the triangular lattice are coupled by near- 
est neighbour 'Heisenberg exchange' coupling that re- 
spects the symmetry of spin-rotations, they order in a 
coplanar 120° ordered state for large spin S, and the 
weight of accumulated evidence [1] suggests that this or- 
dering of Heisenberg antiferromagnets on the triangular 
lattice survives the effects of quantum fluctuations and 
persists down to 5 = 1/2. In the opposite Ising limit, 
when the exchange couplings only couple one component 
of the spins, there are of course no quantum fluctuations; 
however, as was shown long ago by Wannier 0] in the 
5 = 1/2 case, the system remains in a collinear para- 
magnetic state all the way down to T = due to the 
geometric frustration of the triangular lattice. 

Here, our focus is on the role of quantum effects in 
driving the large spin (5 > 3/2) triangular lattice Heisen- 
berg antiferromagnets with relatively strong easy-axis 
single ion anisotropy to an ordered collinear state. As we 
demonstrate below, the leading effects of virtual quan- 
tum transitions induced by the transverse components 
of the exchange coupling pick out an unusual 'orienta- 
tionally ordered' collinear state for all 5 > 3/2 whenever 
the easy axis single ion anisotropy D dominates over the 
isotropic nearest neighbour exchange J. 

In this ordered state which breaks the symmetry of 7t/3 
rotations about a lattice site, but does not break lattice 



translation symmetry (Fig[T]a), the total magnetization 
along the z axis remains zero for small applied field B 
along the easy axis, and the system is thus expected 
to have a m — plateau. This is in sharp contrast to 
5 = 1 magnets with similar easy axis anisotropy — as was 
pointed out in Ref [H, the original results of Refs. [EB0] 
imply that 5=1 antiferromagnets with strong easy axis 
anisotropy on the triangular lattice order in a state that 
has coexisiting three-sublattice spin-density wave order 
for the z component of the spin, as well as spin-nematic 
order in the transverse components of the spins. 

As we detail below, these results are expected to be 
of potential relevance to some members of a class of 
interesting magnets d EH 11, 12, EH in which one 
set of magnetic ions have an unusual trigonal prismatic 
environment that typically leads to significant single- 
ion anisotro py effects in the low energy spin Hamilto- 
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Model and Motivation: We focus on the effects of a 
large easy axis single-ion anisotropy on triangular lattice 
antiferromagnets with Hamiltonian 
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where J > denotes the nearest neighbour antiferromag- 
netic spin exchange interaction between the 5 > 3/2 ions, 
D is the strength of the easy axis single-ion anisotropy 
along the z axis, and B represents the external field along 
the easy axis. 

Because of the layered nature of triangular magnets, 
uniaxial single-ion anisotropy that picks out a common 
easy axis is a generic possibility in such systems. For in- 
his pos sibility is frequently realized on the A 
EE M E3- in the A3ABO3 family of com- 
Ca, Sr, or Ba) 0,0, El El El which 
consist of a triangular array of chains of alternating 'A' 
and 'B' type sites that can host a wide variety of mag- 
netic or non-magnetic cations [Io| with typically antifer- 
romagnetic inter-chain coupling and intra-chain coupling 
of variable sign and magnitude. Although our two dimen- 
sional model does not apply directly in cases with large 



stance, this 
site 0. 
pounds (with A 



2 




FIG. 1: (color online), a) Orientationally ordered state on 
the triangular lattice. The average Ising exchange energy of 
the blue (dark) bonds is lower than that of the green (light) 
bonds, b) The representation of minimally frustrated Ising 
configurations in terms of dimers (dark bonds) on the dual 
honeycomb net. 

intra-chain coupling (such as CaaC^Og), it is neverthe- 
less expected to be a good starting point for other com- 
pounds of this type in which the B site is non-magnetic, 
and the coupling between successive high-spin A sites 
along the chain is concomitantly weak compared to the 
in-plane exchange coupling. 

Returning to our model Hamiltonian, we note that the 
anisotropy term can induce a collinear state for not- very- 
large values of D/J due to the frustrated nature of the 
exchange (for instance, in the classical large 5" limit, the 
system prefers a collinear state for D > 3 J/2). The de- 
tailed nature of the collinear ordering is then controlled 
by the leading effects of quantum fluctuations induced by 
the exchange coupling J. In order to understand these 
effects, we develop a systematic large D/J expansion for 
the leading order effective Hamiltonian that governs the 
effects of these fluctuations — this perturbative analysis is 
expected to yield reliable results whenever the easy axis 
anisotropy forces the spins to be collinear. An analysis 
of the low temperature properties of this effective Hamil- 
tonian then yields detailed predictions for the collinear 
ordering selected in this regime. 

Method: We begin by splitting the Hamiltonian as H = 
Ho+H, with H = -D^XS*) 2 andffx = J£ fe) S v S r 
Clearly, the ground state manifold of Hq at large D is ob- 
tained by requiring all spins to be maximally polarized 
along the z axis and have Sf — o~iS, with o~i = ±1. 
The low energy physics in this Ising subspace is then 
controlled by the leading perturbative effects of the ex- 
change coupling J. These perturbative effects can be rep- 
resented as an effective Hamiltonian H. acting within this 
low energy subspace. Using standard degenerate pertur- 
bation theory and including terms up to 0(J 3 / D 2 ), we 
find (upto unimportant additive constants) for S = 3/2: 

U = Jl ^2 a i a i ~ J 2 J! ~ K^( a i H i + a 3 H i) 

+ J3j2( a ? a ~ +hx -) ( 2 ) 

(ij) 




FIG. 2: (color online), a) Orientational order parameter $ as 
a function of the inverse temperature, b) The double peak in 
the histogram of $ provides a clear signature of a first order 
jump in <E> at the transition. 

where J 1 = ^, J 2 = §jj£, and J 3 = for S = 3/2, 
while the exchange field Hi = V^aj with Fy = 1 for 
nearest neighbours and zero otherwise. 

In the above, the first term corresponds to the leading 
effect of the z component of the spin exchange, while the 
second term arises from the effects of virtual quantum 
transitions of pairs of anti-aligned spins out of the low 
energy Ising subspace: such virtual transitions cost dif- 
fering amounts of energy depending on the local environ- 
ment of the fluctuating pair, and these differences lead to 
different 'zero-point' energies which are encoded in this 
'potential energy' term. Both these diagonal terms gener- 
alize straightforwardly to arbitrary S > 3/2 and we find 
J\{S) = JS 2 and J2(S) = ^F^s-iy 1 • ^ n con trast the 
presence of the off-diagonal term (representing real quan- 
tum transitions) at order 0(J 3 /D 2 ) is special to S = 3/2, 
and it is easy to see that the effects of real quantum tran- 
sitions are in general of order 0{J 2S / D 2S ~ V ) and can be 
safely ignored for S > 3/2. Moreover, even for S = 3/2, 
the J3 term has a significantly smaller numerical coeffi- 
cient compared to the other subleading term J2, and we 
therefore expect the effects of the potential energy term 
J2 to dominate over the off-diagonal J3 term for S = 3/2 
magnets as well. 

Thus, the low temperature behaviour is well-described 
by an effective classical triangular lattice Ising model 
written in terms of Ising pseudospin variables a that 
interact with a leading nearest-neighbour antiferromag- 
netic coupling Ji, and a sub- leading multi-spin interac- 
tion J 2 arising from virtual quantum transitions. At 
low temperature, the leading term constrains the system 
to lie within the macroscopically degenerate manifold of 
minimally frustrated states (in which every triangle has 
precisely one frustrated bond corresponding to a pair of 
ferromagnetically aligned a) of the classical triangular 
lattice Ising antiferromagnet. Our focus here is therefore 
on the possible effects of the subleading multispin inter- 
action term J2 acting within this degenerate subspace. 

Further analysis of the effects of this unfamiliar multi- 
spin interaction on the low temperature physics is facil- 
itated by noting that the minimally frustrated states of 
the classical Ising antiferromagnet may be represented as 
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dimer covers of the dual honeycomb lattice, with a dimer 
placed on every dual lattice link that crosses a frustrated 
bond of the physical lattice (Fig [1] b) . The multispin in- 
teraction J 2 acting within this space of states can then be 
written in this language as an interaction term V between 
the dimers. 

To see this, we first note that this term only acts on 
unfrustrated bonds of a pseudospin configuration. Now, 
since the number of unfrustrated bonds in any minimally 
frustrated configuration is a constant equal to n t , the 
number of triangles in the system, we may rewrite the po- 
tential energy V as: V = —2J 2 (S)n t Y^ n =2 n f™ wnere fn 
denotes the fraction of unfrustrated bonds with n dimers 
on the perimeter of the corresponding double-hexagon of 
the dual honeycomb lattice (Fig[T]b). 

Semiclassics: In order to test this general result for 
the potential energy term V, we have also performed a 
semiclassical expansion in the large S limit and calcu- 
lated the effective Hamiltonian (restricted, as above, to 
the minimally frustrated subspace of the classical Ising 
antifcrromagnet) at leading O(S) order using the proce- 
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dure outlined in Ref 18|. Expanding this large-.!? semi- 



classical result to leading order in the physical small pa- 
rameter J/D, we obtain (apart from unimportant addi- 
tive constants) 



TLsc 



J A S 
32D 2 



(3) 



This result is, at first sight, quite different from the out- 
come of our perturbative analysis. However, it is not 
difficult to see that TLsc is in fact equivalent to the large 
S limit of TL. To see this, we first note that V/2J 2 (S) just 
counts the total number of dimers on the perimeters of 
double hexagons of the dual lattice. On the other hand, 
the semiclassical potential energy TLsc derived above 
can be written as TLsc = ^J2(S)nf l Y^ rl=0 ni 2 g m where 
rih = n t /2 is the total number of hexagons in the dual 
lattice, and g m denotes the fraction of hexagons with m 
dimers on its perimeter. 

Now, each single hexagon with m dimers on its perime- 
ter can be part of 6— m such double hexagons correspond- 
ing to unfrustrated bonds, and these m dimers will thus 
contribute 6 — m times in this sum. As a result, V can 
equally well be written as V = — 2 J 2 {S)nh Y^L=o ~ 
m)g m , Simplifying and taking the large- S limit, we im- 
mediately obtain (apart from an additive constant) the 
alternate single-hexagon form of the potential energy 
that emerged from the semiclassical analysis. 

Low temperature state: Since the average number of 
dimers on the perimeter of a hexagon is 2, it is clear from 
the single hexagon form of V that the minimum potential 
energy is obtained for all configurations with g 2 = 1, i.e. 
for configurations with all hexagons having precisely two 
dimers on their perimeter. To proceed further and un- 
derstand the nature of the low temperature state chosen, 
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FIG. 3: (color online), a) Contour plot of the pseudo-spin 
structure factor S(q) for size L — 24 at 2/3Js(S) = 2.2 (q x 
and q y are projections of the wavevector q along lattice di- 
rections To and Ti (Fig[T])). The Bragg lines with enhanced 
signal are the signature of the orientational order, b) The 
magnetic susceptibility for fields along the easy axis. Note 
the sharp drop in susceptibility when the system enters the 
orientationally ordered state. 



it is necessary to efficiently simulate the corresponding 
interacting dimer model at low temperature. 

In order to do this, we employ a generalization [2(| of 



the procedure of Rcfs 21 



22j and work with the single 
hexagon form of V. Our conclusions from this numerical 
study are readily stated: We find that below a critical 
temperature T c w 1.6772(5), the system orders into a 
orientationally ordered state in which the mean Ising ex- 
change energy on a link of the triangular lattice depends 
on its orientation, but not its position. [We note paren- 
thetically that a very similar ordered state is also realized 
in another interacting dimer model studied recently in an 
entirely different context [Hj].] 

This is seen in the behaviour of the orientational order 
parameter $ = ^ p — £? p e 2p7 ™/ 3 , where B p denotes the av- 
erage of the Ising exchange energy o~iO~j on all links (ij) 
of the p th orientation (p — 0, 1, 2-Fig[T]a) on the triangu- 
lar lattice (Fig [2] a). From the double peak nature of the 
histogram of the order parameter at the transition, we 
see that the transition has a first order character (Fig [2] 
b). 

In such an orientationally ordered state, the Ising pseu- 
dospins are antiferromagnetically arranged in parallel 
rows oriented along one of the three principal directions 
of the triangular lattice. As each such row can be in one 
of two internal states corresponding to the two antifcr- 
romagnetic arrangements of a on that row, the question 
then arises: Are these internal states (that may be repre- 
sented by a 'parity' variable that takes two values) locked 
into some specific ordered pattern, or are they fluctuat- 
ing randomly? To answer this question, we monitor the 
structure factor of the pseudospins a (which corresponds 
to experimentally measurable leading zz component of 
the structure factor of physical spins S) and note that 
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if these internal parity variables had long-range order at 
some wavevector, one would expect to see a correspond- 
ing Bragg peak in the pseudo-spin structure factor, while 
the absence of any observed Bragg peaks would imply 
either a random glassy pattern of parities or rapidly fluc- 
tuating parities. 

As is clear from Fig|3]a, there is no long range order in 
the internal states of the antiferromagnetic rows. Indeed, 
in our numerics, we find that the antiferromagnetic rows 
fluctuate freely between their two allowed internal states. 
However, these fluctuations are only possible because the 
algorithm uses very efficient non-local multi-spin moves 
to rapidly cquilibriate the system with the correct equi- 
librium Gibbs distribution. With realistic single-spin flip 
or spin-exchange dynamics, we find that the antiferro- 
magnetic rows freeze into a random glassy pattern of in- 
ternal parity states as the equilibrium transition temper- 
ature is approached. Such glassy behaviour due to the 
presence of a diverging time scale associated with the 
formation of extended structures has been discussed at 
length in earlier theoretical studies [23(], and our results 



provide another example of such glassy behaviour in a 
disorder-free system. [We also note parenthetically that 
a similar ordered state, but with the row parities 'ferro- 
magnetically' arranged, is realized in the spin-1 triangu- 
lar lattice antiferromagnet AgNi0 2 but the underlying 
mechanism is likely to be quite different [24j.] 

Another interesting aspect of this orientationally or- 
dered state arises from the fact that the total magnetiza- 
tion (M tot ) as well as the uniform magnetic susceptibility 
X to a magnetic field along the z axis are exponentially 
small deep inside the ordered state (see Fig|3]b). As this 
ordered state is stable to small magnetic fields for which 
the potential energy V dominates over the Zeeman en- 
ergy, this implies the presence of low temperature zero 
magnetization plateau that extends for a range of mag- 
netic fields < \B\ < B c ~ J 3 /D 2 . [For \B\ > B c , the 
Zeeman energy gain dominates over the fluctuation in- 
duced potential energy, and the system is expected to or- 
der in the well-known three-sublattice ordered one-third 
magnetization state of the triangular lattice Ising anti- 
ferromagnet.] 

Discussion: As mentioned earlier, our predictions for 
an orientationally ordered state at low temperature are 
expected to be of relevance to other members of the 
Ca3Co203 family of compounds in which the octahedral 
site is non-magnetic, resulting in weak (in comparison to 
the in-plane antiferromagnetic exchange coupling) cou- 
pling between successive high-spin magnetic ions on the 
trigonal prismatic sites along the chain. The orienta- 
tional ordering preserves the translational symmetry of 
the triangular lattice, but breaks the rotational symme- 
try 7r/3 rotations about a lattice site, and is expected to 
give rise to a characteristic set of high intensity lines in 
the zz component of the spin structure factor, which can 
be probed in neutron scattering experiments. Further- 



more, the ordered state is also associated with a zero- 
magnetization plateau in the magnetization response to 
a field applied along the easy axis, while the ordering 
transition is expected to give rise to a sudden drop in the 
magnetic susceptibility \ as we U as a prominent signa- 
ture in the specific heat. In addition, the ordered state 
is also expected to display signatures of slow glassy dy- 
namics of the spins. Given the large variety of magnets 
that crystallize in this structure and the recent interest 
in this family of materials, we hope that our work pro- 
vides some impetus to experimentally identify examples 
of this interesting physics. 
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